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Abstract. We report an analytical expression for the locations of Lindblad resonances induced by a perturbing 
protoplanet, including the effect of disk gravity. Inner, outer and differential torques are found to be enhanced 
compared to situations where a keplerian velocity field for the dynamics of both the disk and the planet is assumed. 
Inward migration is strongly accelerated when the disk gravity is only accounted for in the planet orbital motion. 
The addition of disk self-gravity slows down the planet drift but not enough to stop it. 


1. Introduction 

Interactions between a gaseous disk and embedded proto¬ 
planetary cores could be decisive to understand the dis¬ 
tribution of semi-major axis, eccentricities and masses 
of extra-solar giant planets (Perryman 2000, Schneider 
2004). Exchanges of angular momentum can efficiently 
operate through the excitation of spiral density waves at 
the sites of Lindblad resonances, leading to inward mi¬ 
gration of solid bodies (Goldreich & Tremaine 1979 & 
1980, Ward 1997). For planets with mass less than about 
a Jovian mass, the interaction between the disk and the 
planet is mostly linear (type-I migration), in contrast with 
the non-linear type-II migration where massive planets 
open a gap. In any case, the drift time-scale is much 
shorter than the time required for the completion of a 
(giant) planet. Several works have recently focused on 
means to slow down or stop migration. The presence of 
a toroidal magnetic field (Terquem 2003), tri-dimensional 
effects (Tanaka, Takeuchi & Ward 2002) or corotation 
torques (D’Angelo, Henning & Kley 2002) could act in 
this way. 

From numerical simulations, Nelson & Benz (2003a & 
b) have shown that disk self-gravity can noticeably af¬ 
fect the drift velocity (even for low mass disks). They 
suggested that even very weak changes of the rotation 
curve induced by disk gravity significantly modifiy the 
location of Lindblad resonances, and subsequently the 
total differential torque. Their conclusions are however 
strongly resolution-dependent. In this short communica¬ 
tion, we clarify the influence of disk gravity on type-I mi¬ 
gration by a semi-analytical approach. In particular, we 
determine analytically for the first time the location of 
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Lindblad resonances modified by disk gravity, and com¬ 
pute the corresponding gravitational torques. Analytical 
techniques generally provide reliable diagnostic tools and 
powerful predictions as they implicitely correspond to an 
infinite numerical resolution. 

In Sect. 2, we derive and discuss a general expression 
for the location of Lindblad resonances as functions of 
the disk surface density profile, orbit of the planet, rela¬ 
tive disk mass and edges. In Sect. 3, we successfully check 
this expression in the simple case of radially homogenous 
disks. The effect of the disk mass on Lindblad torques is 
then analyzed. Finally, we consider the case of disks with 
power-law surface density profiles. We conclude in Sect. 4. 

2. Location of Lindblad resonances 

2.1. Background 

A planet embedded in a gaseous disk exerts gravitational 
torques at sites of Lindblad resonances. For low mass bod¬ 
ies, the disk response to the perturbing point mass po¬ 
tential can be considered as linear, so that torques can 
explicitly be determined (Goldreich & Tremaine 1979, 
Artymowicz 1993). The nominal positions Rh{m) of the 
inner (ILRs) and outer (OLRs) Lindblad resonances asso¬ 
ciated with the m-th order Fourier pattern are found from 
the equation D{Ri2) — 0 with 

D = — m^(H — ripi)^, (1) 

where k is the epicyclic frequency defined by = 4D^ -|- 
2RQ.(Kl/dR, n is the fluid angular velocity and flpi is the 
angular velocity of the planet. Generally, Ri^{m) differs 
from the effective locations of Lindblad resonances (here¬ 
after i?*) where the waves become evanescent. Effective 
resonances are found from the equation ZJ*(i?*) = 0 with 
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(Artymowicz 1993, Ward 1997) 


- -O + 
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m^c: 
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where Cg is the sound speed. If the disk aspect ratio H/r = 
?7 is a constant, we then have Cg = and so 




shift 

ILRs 

OLRs 

contribution 

inward migration 

SRi 

> 0 

> 0 

outwards 

slowed down 

5 R 2 

> 0 

< 0 

inwards 

accelerated 

5 R 3 

< 0 

< 0 

inwards 

accelerated 

total 

> 0 

< 0 

inwards 

accelerated 


Table 1. Expected sign of the three individual shifts, and cor¬ 
responding prediction about planetary migration. 


where e = — 1 stands for the ILRs, e = -fl is for the OLRs, 
and / = ^/T _ 1 _ ^ This expression shows that nominal 

and effective resonances coincide only for low m values. 

2.2. EfFects of disk gravity 

The gas in the disk, whatever its mass relative to the cen¬ 
tral mass, is a source of gravity. It modifies the positions 
of the Lindblad resonances in two ways. First, it changes 
the dynamics of both the gas component itself (i.e. disk 
self-gravity) and the planet according to 

= + nl 


the change of the planet angular velocity due to the disk 
mass. 

The “efficiency” of migration can qualitatively be de¬ 
duced from the shifts 6R/R^, of the ILRs and OLRs since, 
in general, the amplitude of gravitational torques de¬ 
pends strongly on the location of the resonances. Further, 
these torques are the largest for intermediate m values 
~ 10 — 20 (~ l/r;). In these conditions, we have / ~ 1, 
F ~ l/me. Using the (crude) monopole approximation for 
Uj (Mineshige & Umemura 1996) and assuming a power 
law surface density prohle, we And 



R* X 




U + 


(4-|-7)e 


—e 


-1 


(9) 


where is the contribution due to the disk, and = 
Uj(i?pi). Second, it shifts the location of effective res¬ 
onances closer to the planet (density waves can propa¬ 
gate between Lindblad resonances). According to Nelson 
& Benz (2003b), the new locations i?dg of effective ILRs 
and OLRs modified by the disk gravity are determined by 
the equation Ildg(Rdg) = 0 where 


Ddg — R* 


2nGTim 

R 


and S is the disk local surface density. 


(5) 


2.3. Second-order differential method 


Provided the relative shift 6R = Rdg — R* is small com¬ 
pared to i?*, we can derive a reliable analytical expression 
for Rdg by a low order expansion of Ddg{R*-\-SR). Keeping 
terms of second order in Ud/Ux only, we find after some 
algebra 


where 


SR = 

SRi + SR 2 + SR 3 

(6) 
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All quantities in Eq. (7) are to be expressed at R*. The 
first two terms in the right-hand side of Eq.( 6 ) are due to 
the disk self-gravity, whereas the third one corresponds to 


Note that I 7 I is of the order of unity in most cases. 
To conclude on the influence of disk gravity, one must 
compare each term in Eq.(9) together. Table 1 summarizes 
the sign of each shift, and the consequence on migration. 
We see that if self-gravity is neglected (i.e. assuming SR « 
SR 3 ), then inward migration is accelerated. This effect has 
been clearly observed in the hydrodynamical simulations 
by Nelson & Benz (2003a & b). On the contrary, self¬ 
gravity alone (i.e. SR « SRi + SR 2 ) slows down migration 
due to an asymetrical shift of the resonances. Further, 
since 


SRi + SR 3 
SR 2 


(5 + 7) 

m 


< 1, 


( 10 ) 


we conclude that the total shift is dominated by the term 
SR 2 . From Eq.(9), we have 


SR 

~R 


OLR 


SR 

~R 


ILR 



4Rpi d 
9m dR 


( 11 ) 

< 0 


For most disks of astrophysical interest, we expect that 
Hd/^K lo be an increasing function of the radius. It then 
follows that OLRs should be shifted more than ILRs, 
meaning that migration tends to proceed inward more 
rapidly when disk gravity in included. 


3. Influence of the disk gravity on low mass 
planet migration 

3.1. Gravitational acceleration in an homogeneous disk 

As Eq.(7) shows, resonance shifts due to the disk mass can 
be determined once the radial gravity held ga = —O^R 
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is known. Since there is no reliable formula for poten¬ 
tial/density pairs in flat disks, we shall consider a simple 
case which allows some analytics, namely a thin disk with 
uniform surface density S (a case often considered in sim¬ 
ulations) . Then, the radial held gji inside the disk exactly 
writes (e.g. Durand 1964) 


gR = 4GS 


E(t’out) I^(t^out) 


^out 


-I- K(Min) - E(Min) 


( 12 ) 


where K and E are the complete elliptic integrals of the 
flrst and second kinds respectively, Uin = a-m/R < 1, 
fout = R/aout < 1, where Oin > 0 is the disk inner edge 
and Oout > flin is the outer edge. A more tractable expres¬ 
sions for gn can be derived from truncated expansions of 
the complete elliptic integrals. For instance, with a classi¬ 
cal second-order expansions over the modulus a; < 1 (e.g. 
Gradshteyn & Ryzhik 1994), namely 


K(a:) = I (1 + + 0(a:4), 

E(a:) = 1(1 - + 0(a:4), 


the held far from edges is given in a good approximation 

by 

gn « -ttGE (vout - wfn) ■ (14) 

Note that contrary to Eq.(12), that formula remains fl- 
nite at the edges, and then is more realistic although ap¬ 
proximated. Since the gravitational potential of the disk 
is minimum very close to the inner edge (like in most as- 
trophysical disks), we can simplify Eq.(14) by considering 
only the term linear in R, and so 

^out 


3.2. Shifts and torques at the Lindblad resonances 


In the homogeneous disk model, we thus have 

‘k “out 


(16) 


where g = M^/M, being the disk mass, and assuming 
<C flout- Hence, 7 = 0. Using Eq.(3), our formula for 
the relative shift becomes 


5R ^^Rl, / m + c/ y A-k;^-2m^ \ 
A* 3flo„t \ m y 1 1 -I- me/ -I- m^r/^ I 


We see that the relative shift is strongly sensitive to the 
location of the planet with respect to the disk outer edge. 
Figure 1 shows as predicted by Eq.(17) for typi¬ 

cal parameters. The agreement between the above formula 
and values determined numerically from Eq.(5) by stan¬ 
dard root finding techniques is excellent. Both Lindblad 
resonances get closer to the planet. The OLRs are slightly 
more shifted than the ILRs, specially for small m (depend¬ 
ing on the disk aspect ratio). This can be understood by 
computing the relative displacement 


5R 

~R 


OLR 


5R 

~R 


ILR 


G _ mflout \ 

V y 


^0. (18) 



Fig. 1. Relative shift SR/R* of the effective ILRs and OLRs 
due to the disk mass, for the following parameters: ain = 0.01, 
flout = 1, p = 0.05 and Rpi/flout = 0.5. 




Fourier pattern m Fourier pattern m 

Fig. 2. Inner and outer Lindblad torques (left) and differential 
torques {right) when the disk mass is accounted for, compared 
to the case without disk gravity. The conditions are the same 
as for Fig. 1. 


Figure 2 displays the Lindblad torques computed following 
the procedure described in Ward (1997). It confirms that 
the outer torques have larger amplitude than inner ones 
when the disk gravity is accounted for. Differential torques 
are larger, and inwards migration should be accelerated. 

3.3. The case of non homogeneous disks 

We have computed the resonance shifts and associated 
Lindblad torques for disks with power law surface den¬ 
sity profiles (i.e. E oc i?“®), typical from disk models and 
observations. The disk gravity field g^ has been deter¬ 
mined numerically from the splitting method described in 
Hure & Pierens (2004). As shown in Fig. 3, resonances 
shifts are weakly affected by the surface density profile 
for s = 1/2 and 1. For steeper profiles however (like for 
s = 3/2), shifts are smaller. Figure 3 also displays the dif- 
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Fourier pattern m Fourier pattern m 

Fig. 3. Relative shifts of the OLRs and ILRs (left) and differ¬ 
ential Lindblad (right) for different power law surface density 
profiles (see text). The conditions are the same as for Fig. 1. 


ferential torques. These conserve the same shape as in the 
homogeneous case, are always larger than in the case with¬ 
out disk gravity. Their magnitude decreases as the surface 
density profile gets steeper, especially for low m values. 
Changes are minor for high Fourier modes. Globally, the 
conclusions established in the homogeneous disk model 
still hold: inward migration should proceed faster due to 
the disk gravity. For large s-values, the effect of disk grav¬ 
ity on migration is predicted to be less and less efficient. 




Fig. 4. Total differential torque versus the planet position rel¬ 
ative to the outer edge for fi — 0.05 (lejf), and versus the 
relative disk mass jj, for Rpi — 0.5 (right), for different power 
law surface density profile in the disk (see text). 

Figure 4 displays the differential torque as functions 
of Rp\ and /t for various exponents s in the power law 
surface density. We see that torques increase as the planet 
orbits closer to the outer edge, and as the disk mass rises. 
Both effects are due to the fact that the the resonances 
get closer to the planet as fj, or i?pi/aout increases, however 
with a slightly larger shift of the OLRs with respect to the 
ILRs (for reasons explained above). 


4. Concluding remarks 

In this paper, we have reported a general expression for 
the shift of the Lindblad resonances due to the disk grav¬ 
ity, whatever the surface density profile, and computed 
associated torques. In contrast with current numerical 
simulations (Nelson & Benz 2003a & b), our analysis is 
not resolution-dependent, thereby enabling reliable pre¬ 
dictions about the migration mechanism of low mass em¬ 
bedded objects. We have considered the effect of the disk 
gravity i) on the planet dynamics, and ii) on the disk itself 
(i.e. self-gravity). Both effects are important and act in op¬ 
posite ways. We confirm that disk gravity plays an impor¬ 
tant role on type-I migration, even for low mass disks. We 
find that the position of the resonances are significantly 
modified and get closer to the planet when the disk mass is 
taken into account. The differential Lindblad torques are 
stronger than in the case where the disk mass is neglected 
(Ward 1997). Our results are also compatible with the 
recent simulations by Nelson & Benz (2003a & b): migra¬ 
tion is accelerated when the disk gravity is accounted for 
in the motion of the planet only, and slowed down when 
self-gravity is added, but does not stop it (assuming that 
all Fourier modes exist). Regarding extrasolar planets, our 
conclusions reinforce the necessity to seek for mecanisms 
able to cancel inward migration. We note that the possible 
suppression of low m-modes of the OLRs (for instance if 
the planet evolves too close to the disk outer edge) could 
decrease the total torque exerted on the planet and change 
its drift. 

It would be of great interest i) to seek for a general 
expression of SR as a function of the surface density profile 
(for instance as an explicit function of the s-exponent), 
and ii) to compare the location of Lindblad resonances and 
associated torques as predicted here with those obtained 
directly by numerical simulations of fully self-gravitating 
disks. This would require very high numerical resolutions. 
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